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We are here interested in i/(z/)-nucleus scattering reactions in the quasielastic (QE) region, where 
I the neutrino interacts with one single nucleon and the mechanisms like RPA or the excitations of 

the target nucleus become less important. In Sec. II a Green's function approach to describe in- 
clusive charged-current (cc) neutrino scattering is developed [1]. In such a process, i^(p)'s interact 
I with nuclei via the exchange of weak-vector bosons and charged leptons are produced in the final 

' state. In our model the conservation of flux is preserved and final state interactions are treated con- 

. sistently with an exclusive reaction. In Sec. Ill a relativistic distorted-wave impulse-approximation 

(RDWIA) calculation of semi-inclusive neutral-current (nc) u- and z^-nucleus reactions is presented 
[2]. The sensitivity to the strange quark content of the nucleon weak current is discussed. 
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INTRODUCTION 



II. THE RELATIVISTIC GREEN'S FUNCTION METHOD 

The cross section of an inclusive reaction where an incident neutrino or antineutrino interacts 



P5 . with a nucleus and only the outgoing lepton is detected is given by the contraction between the 

lepton tensor [1] and the hadron tensor, whose components are given by bilinear products of the 
transition matrix elements of the nuclear weak charged-current operator J'^. The current operator 
is assumed to be adequately described as the sum of single-nucleon currents, corresponding to the 



I weak charged current 

= [fV^m + i^F^^^^-^q^ _ g'at'^t' + Fpq>^j'] r± , (1) 

where are the isospin operators, k is the anomalous part of the magnetic moment, qf^ = {uj, q), 
with = |qp — uj"^, is the four-momentum transfer, and cr^*^ = {i/2) [7^,7^^]. and F2 are the 
isovector Dirac and Pauli nucleon form factors. Ga and Fp are the axial and induced pseudoscalar 
form factors, which are usually parametrized as 

Ga(Q') =5a(1 + QV^a)"' , Fp{Q^) = 2MGA{ml + Qy\ (2) 

where qa = 1.267, m.,^ is the pion mass, and Ma — (1.026 it 0.021) GeV is the axial mass. 
Performing the contraction between the lepton and hadron tensors, the inclusive cross section for 
the QE z^(p)-nucleus scattering can be written as [3] 
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VqRoO + VzzRzz - VqzRqz + VtRt ± V^yRxy , (3) 



where the coefficients v are obtained from the lepton tensor components. All nuclear structure 
information is contained in the response functions R, which are defined in terms of the hadron 
tensor components, i.e., 



Introducing the Green's operator related to the nuclear Hamiltonian iJ, we have 
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for /X = 0, X, y, z, and similar expressions for u^^ and [1]. 

It was shown in Refs. [1, 4] that the nuclear response in Eq. (4) can be written in terms of 
the single particle Green's function, G{E), whose self-energy is the Feshbach's optical potential. A 
biorthogonal expansion of the full particle-hole Green's operator in terms of the eigenfunctions of 
the non-Hermitian optical potential V, and of its Hermitian conjugate V^, is performed 



'S-T- vHe)] I \e)) = , [S-T- V{E)] \ x[ \e)) = . 
Note that E and £ are not necessarily the same. The spectral representation is 
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The hadron tensor components can be written in an expanded form in terms of the single-particle 
wave function, | ipn), of the initial state, corresponding to the energy and whose spectral strength 
is Xn as 
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Tr{£,E) = I 3^\q)^l-V{E) I 

X {x^-\E)\^l-V{E)f{q)\^r.) , 



(8) 
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for /X = 0, X, y, z, and similar expressions for uP^ and u^^ [1]. The factor — V'{E) accounts for 
interference effects between different channels and allows the replacement of the mean field V by 
the phenomenological optical potential Vl- After calculating the limit for rj +0, Eq. 8 reads 



Rerr(-Ef-en,£^f-£n) 
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where V denotes the principal value of the integral. 

Disregarding the square root correction, due to interference effects. The second matrix element 
in Eq. 9, with the inclusion of -v/A^ is the transition amplitude for the single-nucleon knockout 
from a nucleus in the state | ^>q) leaving the residual nucleus in the state | n). The attenuation 
of its strength, mathematically due to the imaginary part of the optical potential, is related to 
the flux lost towards the channels different from n. In the inclusive response this attenuation 



must be compensated by a corresponding gain, due to the flux lost, towards the channel n, by the 
other final states asymptotically originated by the channels different from n. This compensation 
is performed by the first matrix element in the right hand side of Eq. 9, where the imaginary part 
of the potential has the effect of increasing the strength. Similar considerations can be made, on 
the purely mathematical ground, for the integral of Eq. 10, where the amplitudes involved in Tn'^ 
have no evident physical meaning when £ ^ Ef — En- 

In an usual shell-model calculation the cross section is obtained from the sum, over all the 
single-particle shell-model states, of the squared absolute value of the transition matrix elements. 
Therefore, in such a calculation the negative imaginary part of the optical potential produces a 
loss of flux that is inconsistent with the inclusive process. In the Green's function approach, the 
flux is conserved, as the components of the hadron tensor are obtained in terms of the product of 
the two matrix elements in Eq. 9: the loss of flux, produced by the negative imaginary part of the 
optical potential in %, is compensated by the gain of flux, produced in the flrst matrix element by 
the positive imaginary part of the Hermitian conjugate optical potential in x- The cross sections 
and the response functions of the inclusive QE z^(^')-nucleus scattering are calculated from the 
single-particle expression of the hadron tensor in Eq. 10. After the replacement of the mean field 
V{E) by the empirical optical model potential Vi,{E), the matrix elements of the nuclear current 
operator in Eq. 9 are of the same kind as those giving the transition amplitudes of the electron 
induced nucleon knockout reaction [5] and the same RDWIA treatment can be used [6, 7]. 

As a study case, we have considered the ^^O target nucleus and two different values of the 
incident neutrino energy E^, = 500 and 1000 MeV. In order to show up the effect of the optical 
potential on the inclusive reaction, the results obtained in the present approach are compared with 
those given by different approximations. In the simplest one the optical potential is neglected and 
the plane wave approximation is assumed for the final state wave functions x^^^ and x^"-*. In this 
plane wave impulse approximation (PWIA) FSI between the outgoing nucleon and the residual 
nucleus are completely neglected. In another approach the imaginary part of the optical potential 
is neglected and only the real part is included. This approximation conserves the flux, but it is 
inconsistent with the exclusive process, where a complex optical potential must be used. Moreover, 
the use of a real optical potential is unsatisfactory from a theoretical point of view, since the optical 
potential has to be complex owing to the presence of open channels. In Fig. 1, the differential cross 
sections of the ^^0(f^,/Li~) reaction for t?^ = 30 degrees are displayed as a function of the muon 
kinetic energy T^. The behavior of the calculated cross sections is similar for the different energies. 
The effect of the optical potential increases with and decreases increasing E^. The result of the 
PWIA is about 20 30% higher at the peak than the one of the Green's function approach. The 
sum of the exclusive one-nucleon emission cross sections is always much smaller than the complete 
result. The difference indicates the relevance of inelastic channels and is due to the loss of flux 
produced by the absorptive imaginary part of the optical potential. In contrast, the cross sections 
calculated with only the real part of the optical potential are practically the same as the ones 
obtained with the Green's function approach. Although the use of a complex optical potential 
is conceptually important from a theoretical point of view, the negligible differences given by the 
two results mean that the conservation of flux, that is fulflUed in both calculations, is the most 
important condition in the present situation. In contrast, significant differences are obtained with a 
real optical potential in the inclusive electron scattering [4]. The cross sections for the ^^0{i>^, 
reaction are also shown for a comparison. They are always much smaller than the corresponding 
cross sections with an incident neutrino. 
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FIG. 1: The differential cross sections of the ^^0(z/^,/Lt~) reaction for = 500 and 1000 MeV at = 
30 degrees. Solid lines represent the result of the Green's function approach, dotted lines give PWIA, 
long-dashed lines show the result with a real optical potential, and dot-dashed lines the contribution of the 
integrated exclusive reactions with one-nucleon emission. Short dashed lines give the cross sections of the 
^^0(1?^, /U"*") reaction calculated with the Green's function approach. 



III. THE NEUTRAL-CURRENT SEMI-INCLUSIVE SCATTERING 

The differential cross section for the neutral-current z/(i/)-nucleus quasielastic scattering is ob- 
tained from the contraction between the lepton and hadron tensors, as in Ref. [3]. After performing 
an integration over the solid angle of the final nucleon, we have 



da G'e^ . i9 



cos^ 



de dl7 dXN 27r2 2 



VoRoQ + VzzRzz-VQzI^z+VTRT^VxyRxy ^^^^i (H) 

(27r)-^ 



where is the lepton scattering angle and £^n(Pn) t^i^ energy (momentum) of the outgoing nucleon. 
The coefficients v and the responses R are obtained from the lepton and hadron tensor components, 
respectively [2]. 

The transition matrix elements are calculated in the first order perturbation theory and in the 
impulse approximation. Thus, the transition amplitude is assumed to be adequately described as 
the sum of terms similar to those appearing in the electron scattering case [5, 6]. The RDWIA 
treatment is the same as in Refs. [6, 7]. The single-particle current operator related to the weak 
neutral current is 

f = fYj'' + i^P2'y^''Q. - Ga^'^j' + i^pg^' , (12) 



The vector form factors can be expressed in terms of the corresponding electromagnetic form 
factors for protons (F?) and neutrons (F^), plus a possible isoscalar strange-quark contribution 
{FD, i.e., 

^V; Pin) =^{FP_ pny /2 - 2 sin^ ^w^f - Fi /2 , (13) 

where +(— ) stands for proton (neutron) knockout and is the Weinberg angle (sin^ — 0.2313). 
The strange vector form factors are taken as [8] 

^^^^ (l + r)(l + Q7M2)2' "^1^^ ^ (1 + t)(1 + QVM2)2' ^'^^ 

where r = Q'^/{4M^), F|(0) = = -{r'^)/Q, and My = 0.843 GeV. The quantity /n^ is the 

strange magnetic moment and (r^) the squared "strange radius" of the nucleon. The axial form 
factor is expressed as 

GAiQ') = ,o , (15) 

where g\ describes possible strange-quark contributions. 

In order to separate the effects of the strange-quark contribution and of FSI on the cross sections, 
it was suggested in Rcfs. [9] to measure the ratio of proton to neutron yields, as this ratio is expected 
to be less sensitive to distortion effects than the cross sections themselves. Moreover, from the 
experimental point of view the ratio is less sensitive to the uncertainties in the determination of 
the incident neutrino flux. In Fig. 2 the ratio for an incident neutrino is displayed as a function 
of the outgoing nucleon kinetic energy both in RDWIA and RPWIA. The ratio is very sensitive 
to g\ and exhibits a maximum at T/v — 0.6 E^,, which, however, corresponds to values where the 
cross sections are small. The effect is sensibly reduced for g\ = —0.10 with respect to 5^ = —0.19. 
An enhancement of the ratio of ~ 15% is produced by FSI. This result is due both to Coulomb 
distortion and to the different coupling of the optical potential with proton and neutron currents. 
It means that the argument of looking for the strange-quark content in this ratio is strengthened 
by distortion, but the possibility to fix the exact value of the contribution is affected by the 
uncertainties due to FSI. 

Finally, in Fig. 3 we compare our results with the data of the BNL 734 experiment [10]. 
Experimental results were presented in the form of a flux-averaged differential cross section per 
momentum transfer squared Q^. Our results are shown with g\ = —0.19 and without the strange- 
quark contribution. We give also the effect of including the strange vector form factors, Ff and , 
with Ff = -(r2)/6 = 0.53 GcV-^, F|(0) = -0.40 [11], and the dependence given in Eq. 14. 
The strange-quark contribution produces an enhancement of the cross sections, which makes them 
slightly higher than the experimental data. The strange weak magnetic contribution decreases the 
cross section, while the axial and weak electric components give an enhancement. 



[1] A. Meucci, C. Giusti, and F.D. Pacati, Nucl. Phys. A 739, 277 (2004). 

[2] A. Meucci, C. Giusti, and F.D. Pacati, Nucl. Phys. A 744, 307 (2004). 

[3] J.D. Walecka, in Muon Physics, Vol. 11, edited by V.H. Hughes and C.S. Wu (Academic Press, New 
York, 1975), p. 113. 

[4] A. Meucci, F. Capuzzi, C. Giusti, and F.D. Pacati, Phys. Rev. C 67, 054601 (2003). 
[5] S. BofS, et al., Electromagnetic Response of Atomic Nuclei, Oxford Studies in Nuclear Physics, Vol. 20 
(Clarendon, Oxford, 1996); S. Boffi, et al., Phys. Rep. 226, 1 (1993). 



5 



- E, = 500 MeV 




100 200 300 400 500 600 700 800 



[MeV] 

FIG. 2: Ratio of proton to neutron total cross sections of the quasielastic scattering on ^^C as a function 

of the incident neutrino (antincutrino) energy. SoHd and dashed lines are the results in RDWIA with g\ — 
and g% = —0.19. Dot-dashed and dotted lines are the same results but for an incident antineutrino. Long- 
dashed line corresponds to neutrino scattering with g\ = —0.10. 
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FIG. 3: Differential cross sections of the quasielastic scattering, flux-averaged over BNL spectrum [10], 

as a function of the momentum transfer squared. The four upper curves are for incident neutrino and the 
four lower ones for incident antineutrino. Solid lines are the results with no strangeness contribution, dashed 
lines with g\ = —0.19, dot-dashed lines with g\ = —0.19 and F^{0) = —0.40, dotted lines with g\ = —0.19, 
J^I(O) = —0.40 and = —{r'^)/6 = 0.53 GcV~^. Experimental data from Rcf. [10]. The errors correspond 
to statistical and Q'^-dependent systematic errors added in quadrature and do not include Q^-independent 
systematic errors. 



